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ABSTRACT 

There is observational evidence for inside-out growth of giant elliptical galaxies since 
(-h \ z > 2 — 3, which is - in contrast to disk galaxies - not driven by in-situ star formation. 

Many of the ~ 10 systems at high redshift have small sizes ~ lkpc and surface 

brightness profiles with low Sersic indices n. The most likely descendants at z — have, 
on average, grown by a factor of two in mass and a factor of four in size, indicating 
r oc M a with a > 2. They also have surface brightness profiles with n > 5. This 
evolution can be qualitatively explained on the basis of two assumptions: compact 
ellipticals predominantly grow by collisionless minor (mass-ratio 1:10) or intermediate 
(mass-ratio 1:5) 'dry' mergers, and they are embedded in massive dark matter halos 
which support the stripping of merging satellite stars at large radii. We draw these 
conclusions from idealized collisionless mergers spheroidal galaxies - with and without 
dark matter - with mass ratios of 1:1, 1:5, and 1:10. The sizes evolve as r oc M a with 
a < 2 for mass-ratios of 1:1 (and 1:5 without dark matter halos) and , while doubling 
the stellar mass, the Sersic index increases from n ~ 4 to n ~ 5. For minor mergers 
of galaxies embedded in dark matter halos, the sizes grow significantly faster and the 
profile shapes change more rapidly. Surprisingly, already mergers with moderate mass- 
ratios of 1:5, well motivated by recent cosmological simulations, give a ~ 2.3 and after 
only two merger generations (~ 40 per cent added stellar mass) the Sersic index has 
increased to n > 8 (n ~ 5.5 without dark matter), reaching a final value of n = 9.5 
after doubling the stellar mass. This is accompanied by a significant increase of the 
dark matter fraction (from ~ 40 per cent to > 70 per cent) within the stellar half- 
mass radius, driven by the strong size increase probing larger, dark matter dominated 
regions. For equal-mass mergers the effect is much weaker. We conclude that only 
a few intermediate mass-ratio mergers (~ 3 — 5 with initial mass-ratios of 1:5) of 
galaxies embedded in massive dark matter halos can result in the observed concurrent 
inside-out growth and the rapid evolution in profile shapes. This process might explain 
the existence of present day giant ellipticals with sizes, r > ikpc, high Sersic indices, 
n > 5, and a significant amount of dark matter within the half-light radius. Apart 
from negative stellar metallicity gradients such a 'minor' merger scenario also predicts 
significantly lower dark matter fractions for z ~ 2 compact quiescent galaxies and 
their rare present day analogues. 
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1 INTRODUCTION 

Merging is a natural process in hierarchical cosmological 
models and plays a si gnificant role for the assembly of mas- 
sive galaxies (e.g. iKauffmannj 1996). The most massive ob- 
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jects are elliptical galaxies, which are considered to start 
forming their stars at a redshift of z ~ 6 in a dissipa- 
tive environment a nd can rapidly become very massive (~ 



10 n M n ) by z = 2 dKeres et al.ll2005l; Khochfar & Silkl 


2006: 


De Lucia et all l2006t Kriek et all 20061: INaab et all 


2007. 


20091: Joung et al.ll2009l:lDekel et al.ll2009l:lKeres et all 


2009; 


Oser et al. 2010; Feldmann et alj|20ld; iDommsuez Sanchez 



Hilz et al. 



l201ll ; iFeldmann et~al]|201ll ; lOser et al.ll2012l ). Their subse- 
quent evolution is not fully understood yet, as these ellipti- 
cals are observed to be already quiescent at z ~ 2, on average 
4-5 times smaller, and a fa ctor of two less ma s sive than their 
low redshift descenda n ts llDaddi et al. 20051; Truiillol 



Longhetti et all 120071; IToft et all 120071; IZirm et al 



Truiillo et all 120071 ; IZirm et al.ll2007l; iBuitrago et allboOS 



van Dokkum ct al. 2008; Cimatti et al 



2006 



2007 



2008; Franx ct al 



2009 
2010 



Bczanson ct al. 


2009 


Whitaker et al. 


2012 



Dokk um et au feoiO; Saglia 



Onoderall2012l ). This rapid 



evolution eventually proceeds in an inside-out fashion in the 
absence of significant star formation an is accompanied by 
an evolution in the light (and eventually mass ) distribution 
dHopkins et al.ll2009l; Ivan Dokkum et al.ll20ld; lAuger et al ' 



201 lL iTiret et a.1.1 12OIII ; ISzomoru et al.ll2012l ; ISaracco et al 
2012J). This growth process is distinctively different from 
the star formation driven ins i de-out growth of disk galax- 
ies dPrantzos fc Aubertl Il995l; iMatteucci fe Francois! 1 19891: 



Mo et all Il998l; iNaab fc Ostrikeij 120061; iGovernato et all 
20071 ; IWang et all 120111 ; iKauffmann et al.l |2012| ). In ad- 



dition, there is growing observational evidence that 
high redshift ellipticals are more flattened with expo- 
nential like surface brightness distributions (n < 4) 
whereas their potential present day massive descendants 



with 


n > 5 (Toft et al 
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Szomoru et all 20121. 

















Direct observations of massive compact and quiescent 
high redshift systems also provide new constraints on the for- 
mation histories of giant elliptical galaxies. At face value the 
dramatic difference in properties of massive quiescent galax- 
ies at high and low redshift rules out a simple 'monolithic' 
form a tion followed by simple p assive evolution (|Kriek et al.l 
120081 ; Ivan Dokkum et "ail |2008| ) . The observations also rule 
out the formation of ellipticals in a single 'disk merger' event, 
a scenario for massi ve ellipticals which might also suffer 
from other problems (IOstrikeijll980l; [Naab fc Burkertll2003l ; 
iBurkert et al.ll200Sl ; lNaab fc Ostriker&2009l ). Even if the pro- 
genitor disks were gas- r ich leading to a compact remnant 
iRicciardelli et alj|20ld ; IWuvts et ai1l20ld ; iBournaud et all 
l201ll ) those would have to evolve further by a separate pro- 
cess (see e.g. iHopkins et alll2009l ). 

In addition to the well studied stellar components there 
is evidence for the existe nce for dark matter within the half- 
light radii of ellipticals (iGerhard et al.l [20011 ; iThomas et al.l 
|2009| . l201ll ; ICappellaril 1201281 ). Gravitational lensing mea- 
surements of massive galaxies in the local universe predict, 
on average, ~ 30 per cent of the matter within the half- 
light radius of pr esent day massive elliptical galaxies be - 
ing dark matter l|Auger et al.1 12010| ; iBarnabe et alj|201ll ). 
At high redshift, however, the situation is more uncertain 
(|Toft et al.1 120121 ') as dark matter cannot be measured di- 
rectly. The naive expectation is that dark matter is much 
less important due to the dissipative nature of the early 
evolution. Still, assuming its existence and collisionless na- 
ture, dark matter will eventually affect the distribution of 



stars during the further assembly of the galaxies (see e.g. 
lLackner fc Ostrikerll201fj| ; Iffilz et al. 1 120121 ) . 

The underlying idealized assumption for the study pre- 
sented here is that present day massive elliptical galax- 
ies (> 10 1 Mq) since z ~ 2 have grown predominantly 
by accreting stars that have formed in other galaxies and 
this process can be approximated by simulations of merg- 
ers of collisionless 'dry' stellar systems. This assumption 
is not too far fetched as there is clear observational evi- 
dence for the predominance of old (z > 2) stellar popu- 
lations in these galaxies, leaving little r oom for gas accre- 



Brinchmann & Ellis 2000; Treu et al.l 120051; IThomas et all 


20ld; 


van Dokkum ct al. 


20ld; Lopez-Saniuan et al.l l20ld; 


Young 


l201ll; ICrockerl 


2012). Directly observed merg- 



ers, inferred merger rates, as well as stellar fine struc- 
ture at large galactocentric radii are indicative of stel- 
lar merger and accretion events dMalin fc Carted Il983l; 
Schweizer fc Seitzerj [l99p: iTran et al l 120051; Ivan Dokkuml 



20051; iBell et al.l lOOGallbl: iFaberl 20071; iMcIntosh et al.HoolI 
Whitaker fc van Dokkuml 120081; iRobaina et al.1 |2010|; iDucl 
. 2011] ; [Lotz et al.l 1201 ll; ITruiillo et al. l 1201 ll ; IWeinzirl et~aH 
20 111 ; (Tal et al.1 12012| ; iLopez-Saniuanl \20m . On the theo- 



retical side, phenomenological as well as direct numerical 
studies predict that the late assembly of massive galaxies is 



(Kauffmannl 19961; 


De Lucia et al.1 20061; 


Ciotti et al.l 


2007 


Naab et al.1 120071; 


Guo & White! 120081; 


Naab et al.l 


20091 


Yang et al.l|201l|; lLackner fc Ostrikerll2010l;IOser et al.1 


20iq 



Guo et al.l 120111; IFeldmann et all 120111; lOser et al.l 1201 
Yang et al.l |2012|; iGabor fc Dave! [20121 ; iMoster et all [201 
Lackner et al.l 20121 )). In a previous study Hilz et al. (|201 



used idealized merger simulations of spheroidal galaxies to 
re-investigate the effect of the galaxy merger mass-ratio - 
limited to the extreme cases of 1:1 and 1:10 - on the merger 
driven evolution of compact high-redshift spheroids. The 
main finding of this study was that both minor and ma- 
jor mergers lead to size growth and an increase of the dark 
matter fraction, however by different physical processes and 
at significantly different strengths. Violent relaxation in ma- 
jor mergers mixes dark matter into the central regions and 
escaping particles limit the expected size growth even be- 
low the expected values. In 1:10 mergers on the other hand, 
satellite particles are stripped at large radii where the initial 
host galaxies are dominated by dark matter. As a result the 
stellar effective radii and the dark matter fractions of the 
galaxies grow rapidly. 

Here we extend the lHilz et al.l (|2012l ) study to focus on 
the inside-out growth including a moderate mass-ratio of 
1:5. This step is important for a number of reasons: The 
average stellar mass growth since z ~ 2 is a factor of two. 
Mergers with mass-ratios of 1:10 typically take a long time 
to complete and it is not clear whether ten of these merg- 
ers can even be completed in a Hubble time. Additionally, 
recent observationally estimated merger rates indicate more 
minor mergers for massive galaxies and at higher redshift 
but the rates are still low and it is not clear - also theoret- 



the size growth dNipoti et al. 009b|la!: Williams et al. 


201 ll; 


Weinzirl et al.ll201ll;lNewman et al.l2012l;lNipoti et al. 
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Cimatti et al.1 |2012|; iMan et al. 


120121; lEdwards & Pattonl 


20121; iMarmol-Queralto et al.l 
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Figure 1. Stellar mass versus effective radius for the initial 
galaxy models. For our fiducial set of simulations (diffuse satel- 
lites) all galaxies have an effective radius of Ikpc (filled cir- 
cles). For the second set of simulations with compact satellites 
(op en circles) the sizes follow the z = 2 relation (red line) 
by IWilliams et al. I j201dl V The initial conditions are well sep- 
arated from local mass-size relations of ellipticals lShen_et_aL| 
20031; iHvde fc Bernardill2009l ; iGuo et al.ll2009l ; iNipoti et alj|009al ; 
Auger et al.ll2010ri . 



mergers result in a more rapid mass growth and fewer of 
them would be needed provided the size grows rapid enough. 
Another reason to focus on mass-ratios of 1:5 was recently- 
provided by direct cosmological simulations. Three indepen- 
dent numerical studies have found that the average mass- 
weighted mass ratio of galaxy mergers building massive el- 
liptical s is 1:4 - 1:5, making this regime particularly inter- 
estin g (|Oser et al.ll2012l : lGabor fc Davell2012l :l Lackner et aJj 
I2012T I. 

In section [2] we give a short review of the initial galaxy 
models and the simulation parameters. In the subsequent 
Sections [3] [4] and [6] we investigate the evolution of sizes, 
surface densities, profiles shapes and dark matter fractions, 
respectively. We discuss the results and conclude in Section 

m 



2 SIMULATIONS 

Continuin g the long li ne o f research in this di- 
rection (Iwhitej Il978l. Il979l ; iMiller fc Smlt 



1980 



Villumsenl Il983l; iFarouki et all 119851: INipoti et al.l 12001 
Bovlan-Kolchin et al.ll2005l ; ICiotti et al.ll2007l ; INipoti et ail 
009b) we present a set of simulations of dissipationless 



mergers of spheroidal galaxies with and without dark mat- 
ter halos and with mass ratios of 1:1, 1:5, and 10:1. Details 
about the initial models and the simul ations parameter s 
for 1:1 and 1:10 mergers are presented in lrlilz et al.l |2012i ). 
The simulatio ns were performed w it h the N-body/SPH 
code VINE (|Wetzstein et all 120091 ; iNelson et all |2009|) . 
For the present study we have also performed comparison 



run with GADGET (|Springelll2005h with identical results. 
Here we only give a brief summary of the simulation 
setup. The initial galaxy models are isotropic, spherically 
symmetric one- and two- component syst ems following 
Hernquist density profiles (|Hernquistl ll99Cf ) either for a 
model representing only the stellar component of the galaxy 
(one-component, bulge-only) or a bulge embedded in a 
massive dark matter halo (two-component, bulge+halo). 
The host bulges have a stellar mass of M t ^ OBt — 1 and a 
scale radius of a^host = 1-0 (correspondin g to a projected 
half- mass radius of R c ~ 1.8a,. host (|Hernauistl Il990l ) 
realized with 100.000 particles. For the two-component 
cases we assume an additional dark matter halo with a 
total dark to stellar mass ratio of A4dm/M» = 10 and a 
ratio of the scale radii of ad m /a* = 11- The halo is realized 
with 1.000.000 particles resulting in equal-mass particles 
for both components. The gravitational softening length is 
set to e = 0.02 for all particles. The stability of the initial 
N-body models was demonstrated in lrlilz et al.l (|2012l ). The 
merging satellite galaxies are 5 (mass-ratio 1:5) or 10 times 
(mass-ratio 1:10) less massive than the host and have a 
correspondingly lower particle number. For our fiducial set 
of simulations we assume the same scale radius, a* jSat = 1.0, 
for all bulges, resulting in satellites being more diffuse than 
the host. For a second set of simulations we chose more 
compact satellites with scale radii of a* jSa t = 0.8 for the 1:5 
mergers and a*, sa t = 0. 5 for the 1:10 merge rs following a 
z ~ 2 mass-size relation (Wil liams et al.|[2010t) . 

The open and filled circles in Fig.[T]indicate the location 
of the initial galaxy models in the mass-size plane compared 
to several published present day and high redshift mass-size 
relations. The initial massive host galaxies resemble typi- 
cal massive and compact 'red nuggets' at z = 2 which are 
smaller t han typical present d ay galaxies at 10 1 Mq (red 
line from IWilliams et al.l 120101 ). However, for redshift two 
galaxies below ~ 4.0 x 1O 1O M0 (which would be the mi- 
nor merger partners) there is no consistent size information 
available. Therefore we assume two scenarios: a fixed size 
of ~ lkpc for the satellite s, similar to local galaxies (see 
e.g. iMiseeld fc Hilkenl201ll ). indicated by the filled red and 
green circles. Here the merging lower mass galaxies are dif- 
fuse spheroids with low phase space densities comparable 
to low mass disk-like galaxies. Alternatively we use smaller 
sizes consistent with a simple continuation of the z = 2 
mass-size relation to lower masses (open circles), resulting 
in compact spheroids. 

The first generation of equal-mass mergers are parabolic 
mergers of one- and two-component models of the initial 
host galaxies. The second generation is a re-merger of the 
duplicated, randomly oriented, first generation merger rem- 
nant, which was allowed to dynamically relax at the cen- 
ter. The randomly oriented galaxies approach each other on 
parabolic orbits with a pericenter distance of half the spheri- 
cal half-mass radius of the progenitor remnants, i.e. the peri- 
center distances increase with each merger generation. The 
sequences of intermediate and minor mergers with initial 
mass-ratios of 1:5 (1:10) are also simulated with one- and 
two-component models. In the following we use the tradi- 
tional convention to call all mergers with mass-ratios smaller 
than 1:4 minor mergers. Initially, the mass-ratio is 1:5 (1:10) 
and the galaxies are set on parabolic orbits. The randomly 
oriented merger remnants of the first generations are then 
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Figure 2. The projected spherical half- mass radius of the stellar 
component (the mean value along the three principal axes) as a 
function of bound stellar mass for 1:1 (blue), 1:5 (red), and 1:10 
(green) mergers . The observed size g r owth is indicated by the solid 
black line (see [van Dokk um et all ||2010| ) for the observational 
uncertainties). The size evolution of models in the grey shaded 
area is too weak to be consistent with observations. All mergers 
of bulges embedded in massive dark matter halos and high mass- 
ratios (1:5, 1:10, red and green solid/dashed-dotted lines) show a 
rapid size evolution. The size evolution of the bulge-only models 
(short and long dashed lines) are not efficient enough, except 
the 1:10 scenario with a diffuse satellite (green dashed line). The 
accretion of compact satellites results in less size growth compared 
to the diffuse satellites. The major merger lines (blue) show size 
growth that is far too slow. 



set on parabolic orbits with the initial satellite galaxy mod- 
els and a mass-ratio of now 1:6 (1:11), and so on. We per- 
form 6 generations of 1:10 mergers and 5 generations of 1:5 
mergers using satellites with fixed scale radii (Sat 1:5/1:10, 
see Fig. [TJ. For comparison, we also performed minor merg- 
ers of two-component models with compact satellites (cSat 
1:5/1:10, see Fig. [I}. Again all parabolic orbits have peri- 
center distances of half the spherical half-mass radius of the 
bulge of the massive progenitor galaxy. 

The 1:10 bulge-only simulations have no dark matter 
halo and the in-falling satellites suffer less from dynami- 
cal friction. Therefore the final coalescence takes by far the 
longest time. However, taking about ~ 9Gyrs for 10 merger 
generations even this process could be completed by the 
present day assuming a z = 2 progenitor. All other merger 
series are completed in less than ~ 7Gyrs. 



size p lane from z ~ 2 to the present day (|van Dokkum et al.l 
2010). The models here are idealized in the sense that they 
do not contain any dissipative component. If present, like 
in the 'real' universe, thi s component would reduce the size 
growth per added mass jPekel fc Coxi r2006). Therefore we 
consider a model a failure if it occupies the the shaded area 
in Fig. [T] Promising models have to lie above this line so 
that small amounts of gas - and therefore somewhat smaller 
sizes - can be tolerated. 

Equal-mass mergers show an almost linear increase 
of size with mass, (see also Bov lan-Kolchin et al.l 12005 



Ciotti et al1l2007l: iBezanson et alj 2009; Ni poti et aLll009b 



Hilz et al. [2012). independent of whether the stellar system 



is embedded in a dark matte r halo or not (blue so l id an d 
dashed lines). As discussed bv lBovlan-Kolchin et al.l (|2005l ). 
in mergers with dark matter halos the in-falling galaxy suf- 
fers more from dynamical friction in the massive dark mat- 
ter halo of the companion galaxy, resulting in more energy 
transfer from the bulge to the halo, leading to a more tightly 
bound bulge with a smaller size (blue solid line, Fig. [2j com- 
pared to the model without dark matter (blue dashed line, 
Fig- HJ. If we combine the results of both major merger sce- 
narios this yields a mass-size relation of r e ocM 0,91 whic h 
is comparable to the results of lBovlan-Kolchin et al.l l|2005h , 
who found a slightly smaller exponent (~ 0.7) for orbits with 
high angular momentum and an exponent ~ 1 for pure ra- 
dial orbits. Nevertheless, as the size grows at most linearly 
with mass, dissipationless major mergers cannot be the main 

driver for the size evolution of early-type galaxies. 

As expected from s imple virial estimates (ICole et al 



l200d: ICiotti et all 120071 : iNaab et all 120091: IBezanson eta! 



120091 ). the size evolution is stronger for bulge-only models 
with lower mass-ratios of 1:5 and 1:10 (red and green lines 
in Fig- [2J - However, except for the 1:10 mergers with a dif- 
fuse satellite (green short dashed line), all minor mergers 
with bulge-only satellites are not efficient enough to escape 
the 'forbidden' area. This behaviour is improved for minor 
mergers of two-component models, where bulges are embed- 
ded in a massive dark matter halos. For mass-ratios 1:5 and 
1:10 the size evolution is in excess of the observed evolu- 
tion. In the case of 1:5 minor mergers with a less compact 
satellite (red solid line), we obtain a mass-size growth rela- 
tion of r e ocM 2,4 with a similar or even larger exponent for 
both two-component 1:10 scenarios. Therefore we consider 
all minor merger models (1:5 and 1:10) with dark matter 
halos and the diffuse 1:10 mergers to be consistent with ob- 
servations even in more realistic mod els, where dissipational 
effects would reduce the size growt h (Robertson ct al. l2006l : 
ICox et al.ll2006l : iHopkins et aLll200Sl : ICovington et al.ll201ll ) . 
The size growth via major mergers is too slow to fit obser- 
vational data. 



3 EVOLUTION OF THE SIZES 

After the completion of every merger, we allow the central 
region of the remnant to relax, before we compute the pro- 
jected circular half- mass radii, r e , along the three principal 
axes and the bound stellar (bulge) mass, M*. 

In Fig. [2] we show the evolution of the half-mass radius 
as a function of the bound stellar mass for 1:1 (blue), 1:5 
(red), and 1:10 (green) merger hierarchies. The black line 
indicates the observed evolution, r e ocM 2 ' 04 , in the mass- 



4 EVOLUTION OF THE SURFACE DENSITIES 

We now take a closer look at the evolution of the surface 
densities of the merger remnants. In Fig. we plot the stel- 
lar surface mass densities and cumulative stellar mass pro- 
files after every 1:1 and 1:5 merger and every second 1:10 
merger for bulge-only models (top panels) and bulge+halo 
models (bottom panels). For equal- mass mergers there is 
mass growth at all radii, i.e. the lines are shifted more or 
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Figure 3. Top panels: Stellar surface mass density profiles and cumulative stellar mass distributions for the fiducial 1:1 (left), 1:5 
(middle) and 1:10 (right) bulge-only merger models after each merger event (colored lines from bottom to top, for 1:1 and 1:5 mergers 
we show every generation, for 1:10 mergers every second). For all mass-ratios the surface densities increase at all radii with merger 
generation with a trend of a stronger increase in the outer regions for lower mass-ratios. Bottom panels: Same for mergers of systems 
with dark matter halos. Equal-mass mergers show a similar behaviour than bulge-only models. For 1:5 and 1:10 bulge+halo models 
(mergers with compact satellites are indicated by the dotted lines) the inside-out growth trend is stronger than for bulge-only models. 
Regions inside 4kp c are a lmost unaffected and the satellite stellar mass assembles predominantly at large radii similar to observations 
Ivan Dokkum et alj <201C|) . 



less parallel to high er densities independent of the presence 
of a halo (see e.g. iMiller fc Smith! Il980l : IVillumsenl Il983l: 
iFarouki et al, I ll983l for discussions on the weak but exist- 
ing break of homology) . This ev olution scenario wou l d be i n 
disagreement to observations of Ivan Dokkum et"aH (|2010l ). 
which show, that early-type galaxies grow inside-out, i.e. 
the central densities stay constant and most of the mass as- 
sembles at larger radii, building up an extended envelope of 
stars. 

The second column in Fig. [3] depicts the surface densi- 
ties and mass assembly of minor mergers with an initial mass 
ratio of 1:5. For the bulge-only models (top) with a diffuse 
satellite (Sat 1:5, Fig. [2}, the surface density stays nearly 
constant up to r ~ 1 — 2kpc and increases mainly in the 



outer parts. This effect was already reported by IVillumsenl 
( 1983) as a possible explanation for abundance gradients 
in elliptical galaxies. The same scenario leads to an even 
stronger inside-out growth if the galaxies are surrounded by 
a dark matter halo (bottom panel, second column) . Here the 
bulge particles get stripped at larger radii and the central 
surface density (r < 4fcpc) stays unaffected; it only increases 
at radii r > 2 — 3kpc. The size growth shown in Fig. [2] is 
due to this build-up of a massive stellar envelope. The mass 
is added to large radii where, prior to the merger, the dark 
matter component was large. Due to our accounting pro- 
cedure the dark matter fraction within r c increases, even 
though the dark matter added in such mergers to the in- 
ner parts is negligible. The dotted lines in these panels show 
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Figure 4. Surface brightness profiles fiy(r) of the fiducial two-component 1:1 (left), 1:5 (middle), and 1:10 (right) merger remnants as a 
function of radius for the initial conditions (black circles), the first merger generation (blue circles, for 1:10 we show the second generation) 
and the final remnants (red circles). The over-plotted dashed lines show the best fitting Sersic function for the outer profiles. The data 
was fitted from 0.02r e to either 10r e or a limiting surface brightness of my = 27mag/arcsec 2 . The residuals Afi < 0.2mag/arcsec 2 (lower 
panels) are small, except for the 1:10 case, which cannot be reasonably well fitted anymore for late generations (n > 20) The fitted 
effective radii r G g t (short vertical lines) are slightly smaller than the projected half-mass radii r c (arrows). For minor mergers - 1:5 and 
1:10 - the best fitting Sersic index n increases rapidly as the mass added in minor mergers grows. Note that for minor mergers the central 
surface densities decline the most even as the outer surface densities (and Sersic indices) increase. 



the four remnants, where the satellites are more compact 
(cSat 1:5, Fig. [2]) and lie an extension of the z ~ 2 mass-size 
relation. Obviously, the results are very similar. 

The six generations of minor mergers with an initial 
mass-ratio of 1:10 are shown in the right column of Fig. [3] In 
the case of bulge-only models (top panels) , the surface den- 
sity increases predominantly at larger radii, similar to the 
previous scenario, but now the satellite stars are even less 
bound compared to the 1:5 case and therefore are stripped 
at larger radii, even without a dark matter halo. If present, 
this effect is enhanced (lower panels of right column), as the 
satellites first orbit through the massive dark matter halo - 
and lose stars - before they reach the center of the host. The 
surface densities are unaffected inside r = 5kpc. In sum- 
mary, 1:5 and 1:10 mergers lead to a significant change in 
the mass distribution of the galaxies with most of the stel- 
lar satellite material assembling at large radii. This picture 
hardly changes for compact satellites (cSat, dotted lines), 
although the scale lengths of the satellite stars are signif- 
icantly smaller, they are more bound and resists the drag 
force of the host potential for a longer time. Consequently, 
somewhat more material gets closer to the central regions. 
In summary a minor merger driven evolution scenario - in 
particular for 1:5 bulge+halo mergers - is in good agreement 
with the picture we get from observations. Although we are 
definitely affected by poor numerical resolution at the very 
centers of the galaxies we note the differential effect that 
for minor mergers the central surface densities decline the 



most even as the outer surface densities (and Sersic indices) 
increase in agreement with expectations from rec ent obser- 
vations (|Szomoru et al.ll2012l : iTruiillo et al.ll2012T ). 



5 EVOLUTION OF PROFILE SHAPES 

The curvature of the light profiles of elliptical galaxies is 
an important parameter as it correlates with other ob- 
served properties of elliptical galaxies, such as the effec- 
tive radius r e , the total luminosity and the stellar mass 
ilCaon et al.lll993l: iNipoti et alj|2003 : iNaab fc Truiilloll2006l : 
iKormendv et al] 120091 ). We therefore use a Sersic r 1 ^" 
( Sersidll968l ) function to fit synthetic surface brightness pro- 
files of our simulations, 

J(r) = I < ,.10- 6 »« r / r -> 1/n - 1 \ (1) 

where the three free parameters are the effective surface 
brightness I e , the effective radius r e and the so called Ser- 
sic index n. Here, n = 1 corresponds to an exponential 
profile and n = 4 to the familiar de Vaucouleurs profile 
(|de Vaucouleursll94Sl ) . The factor b„ , is chosen such that the 
effective radius r e encloses half of the total luminosity. For 
the expected range of Sersic indices, this factor can be ap- 
proxi mated by the relation b n = 0.868n — 0.142 (|Caon et al.l 
Il99&t ). We convert the projected surface mass densities dis- 
cussed in section [3] to a V-band surface brightness profile as- 
suming a stellar age of 10 10 yr and solar metallicity Z = 0.02 
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l|Bruzual fc Charlotll2003l) which is then fitted with a Sersic 
function (see lNaab fc Truiillol l|2006h for details of the fitting 
procedure). Those population properties are reason able for 
present day massive ellipticals (|Thomas et al.ll2010l ) but not 
necessarily so for higher redshifts. However, details of the 
choice of population age and abundance are unimportant for 
the determination of the profile shape (which are the same 
as for the surface density profiles) as long as we assume a 
constant mass-to- light ratio (for the stellar component) with 
radius. 

Figure [3] shows examples of the Sersic fits to the sur- 
face brightness profiles and the residuals of our fiducial 
bulge+halo mergers with mass-ratios of 1:1 (left), 1:5 (mid- 
dle) and 1:10 (right). The profiles are fitted from 0.027V to 
either 10r e or to a limiting surface brightness of my = 2 7 
mag/arcsec -2 (|Truiillo et al.ll2004l ; IKormendv et al"1l2009l ). 
The residuals are very small (A/j, < 0.2mag/arcsec 2 ), except 
for the innermost regions, where the profiles are affected by 
insufficient sampling, softening and relaxation effects. The 
initial Hernquist spheres (black cir cles in all panels) hav e 
a Sersic index of n = 3.9 (see also iNaab fc Truiillol 12006). 
In all cases the fitted Sersic profiles overestimate the cen- 
tral surface brightness leading to a slightly smaller fitted 
effective radius r Ci g t (narrow vertical lines at the bottom of 
each surface brightness panel) than the projected half-mass 
radius r e (corresponding arrows) presented in Fig. [2] The 
differences are small and do not affect any conclusions in 
the paper. 

In the case of 1:1 mergers of two-component models 
(left panel, Fig[4|, we can see that the profile shape hardly 
changes for the remnants. The Sersic index increases weakly 
from n = 3.9 to n = 5.5 with a corresponding increase of 
mass by a factor of about two, hardly enough to explain the 
very high numbers obse rvers find for large massive ellipti- 
cal g alaxies (n ~ 10, see lCaon et al.lfl993l ; IKormendv et al.l 
120091 ). This picture changes dramatically for minor mergers 
where the Sersic index increases rapidly to values of n ~ 10 
for a similar increase in mass. This is a direct consequence 
of the different evolution of the surface density profiles dis- 
cussed in section [4] 



In Fig. [5] we show, for all fiducial bulge-only and 
bulge+halo models, the evolution of the Sersic index as a 
function of the assembled bound stellar mass. There is only 
a moderate increase in Sersic index for the equal-mass merg- 
ers. In contrast, after two generations of 1:5 mergers with 
a mass increase of of only 40 per cent the Sersic index can 
be n > 7 (for the bulge+halo model) and the final remnant 
reaches values of n ~ 9.5, which is in the rang e of the ob- 
served present day massive elliptical galaxies (|Caon et al.l 
1 19931 : IKormendv et al.ll2009l ). The corresponding bulge-only 
minor merger scenarios (red and green dashed lines in Fig. 
[5]) show similar but weaker trends yielding final Sersic in- 
dices of n ~ 7.5. Still, the overall evolution is much faster 
for two-component models. This again indicates that the 
merger mass-ratio and dark matter halos play an important 
role as they increase the effect of stripping at large radii in 
a way, that the accreted stellar mass assembles at the 'right' 
regions of the host galaxy. 
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Figure 5. Evolution of the Sersic indices for all fiducial ma- 
jor and minor merger remnants shown in Fig. [4] as a function 
of bound stellar mass (the evolution for compact satellites is sim- 
ilar). All equal-mass mergers show a very weak evolution (black 
lines). Lower mass-ratios lead to a stronger evolution for 1:5 (red 
dashed) and 1:10 (green dashed) bulge-only models. The presence 
of a dark matter halos even enhances this effect as more stars are 
deposited at large radii (red and green solid lines). The 1:5 two- 
component model shows the strongest evolution reaching values 
of n > 9 after mass increase of only 50 per cent. For 1:10 merg- 
ers we only show the first generations with reasonable Sersic fits. 
Later generations cannot be accurately fitted anymore (see Fig. 



0. 



6 EVOLUTION OF DARK MATTER 
FRACTIONS 

In this section we investigate the dark matter fractions 
within the effective radii of our simulated merger remnants 
in the light of recent lensing observations, which predict 
an inc reasing dark matte r fraction for more m assive ellip- 
ticals dAuger et al.l l2010l: iBarnabe et al.ll201ll). d ynamical 



modeling (|Gerhard et aki |200ll ; ICappellaril l2012al ) and the 



pos sibly low dark m atter fractions in high redshift galax- 
ies |Toft et al.| [2012). The dark matter fractions f d m for all 
bulge+halo simulations 



fdm(r < r 50 ) = M dm {r < r 50 )/M t ot(r < r 50 ) 



(2) 



are shown in Fig. [6] as a function of assembled stellar mass. 
Here rso denotes the spherical half-mass radius of the stel- 
lar component and Mdm, Mtot are the halo mass and total 
mass within rso. We have already explained how the ob- 
servable /d m will increase simply by accounting, when the 
stars are added far out in the dark matter halo. The dark 
matter fraction increases rapidly with each subsequent mi- 
nor merger generation regardless of the mass-ratio. For a 
mass increase of a factor of two the dark matter fractions 
increases by almost a factor two from ~ 40 per cent to 
> 70 per cent for minor mergers whereas equal-mass merg- 
ers only an increase to ~ 55 per cent. The strong evolution 
with added mass for minor mergers is a consequence of the 
rapid size growth ( Fig. [H, which is in good agreement with 
l|Nipoti et al.|[009b1 ). The evolution of fdm. correlates with 
the radii of the merger remnants. Therefore the 1:5 scenario 
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Figure 6. Evolution of the dark matter fraction, /j m (r < r^o) = 
-Mdm( r < '*5o)/-Mtot(r < ?*5o)) within the spherical half-mass 
radius r'50 (which is similar to the half mass radius) for all 
bulgc+halo models as a function of bound stellar mass of the 
merger remnants. Equal-mass mergers show a weak increase of 
/dm- A very rapid increase of f^ m is found for 1:5 (red) and 1:10 
(green) mergers, with a similar evolution for compact satellites 
(dashed lines). The evolution of /<j m for 1:5 mergers with mass is 
very similar to 1:10 mergers indicating that less events results in 
similar change of /d m . 



with more compact satellites (red dashed line), which shows 
slightly weaker size growth than the fiducial case (Fig. [5}, 
also has a slightly lower dark matter fraction. On the other 
hand, the 1:10 bulge+halo remnants grow rapidly in size and 
therefore have the highest dark matter fractions. It is easy 
to understand this trend. As noted, the dark matter is not 
pushed inwards; rather we are adding stars to the outer, dark 
matter dominated parts of galaxies so naturally the amount 
of dark matter within the stars increases. A detai led analysis 
of this process can be found in Iffllz et all l|2012l) . 



7 DISCUSSION AND CONCLUSION 

We present idealized binary merger simulations of spheroidal 
galaxies with initial mass-ratios of 1:1, 1:5, and 1:10 repre- 
sented by pure bulges and bulges embedded in dark matter 
halos. The 1:5 mass-ratio is of particular interest as it is 
similar to the average mass weighted galaxy merger mass- 
ratio driving the assembly of massive galaxies recently pub- 
lished cosmological simulations and therefore those minor 
mergers (or intermediate mass-ratio mergers, traditionally 
all mergers with mass ratios smaller than 1:4 ar e called mi- 
nor mergers) represent a typical growth mode l|Oser et al.l 
120121 ; iGabor fc Pavel I2OI2I ; lLackner et al.l l2012h . With this 
study we aim at assessing the effect of pure collisionless ma- 
jor and minor mergers on the structural evolution for an 
already existent, observed, high-red shift population of com - 
pact and massive galaxies (see e.g. iBournaud et alj (|2007l ) 
for a discussion on the effect of repeated minor mergers on 
disks). 

Confirming earlier studies we find that major merg- 



ers alone cannot explain the observed size growth but 
minor mergers with mass-ratios of 1:5 or 1:10 show a 
significantly stronger size growth per added stellar mass 
llKhochfar fc Silkll2006l; INaab et aLll2009l; iNipoti et alJl009al : 
IHopkins et al.ll2010l ; iHilz et al.ll2012l : lOser et al.ll2012l ). This 
effect is slightly enhanced for more diffuse satellites and sig- 
nificantly enhanced if the galaxies are surrounded by dark 
matter halos. In these cases the satellite stars are more effi- 
ciently stripped at larger radii, either because they are less 
bound or because they orbit in the deeper pote ntial wells 
of the dark matter halos. The latter effect (see IHilz et all 
l|2012h for a detailed discussion) leads to an inside-out size 
growth of r oc M 2 ' 3 already for merger s with a mass-ratio 
of 1:5 , in agreement with o bservations (|van Dokkum et al.l 
l2010l : ISzomoru et all 120121 ). Overall, the picture is very 
similar to the dy namical friction driven gala ctic cannibal- 
ism described by lHausman fc Ostrikerl (Il978l) as an expla- 
nation for BCG p roperties l|Ruszkowski fc Spri ngcl 2009; 
lLaporte et~ai1l2012l ). 

The inside-out assembly of mass at larger radii results 
in a significant change of the surface density profiles which is 
here quantified by an increase of the Sersic index n. Starting 
at a fiducial value of n = 4 the increase is weak for major 
mergers but very strong for minor mergers of galaxies em- 
bedded in dark matter halos. In the most extreme case only 
two 1:5 mergers of bulges with dark matter halos change the 
Sersic index from n ~ 4 to n ~ 8.5 at a concurrent stellar 
mass increase of only 40 per cent. We note that n = 4 is the 
fiducial starting point of the specific simulation setup we 
have chosen. Qualitatively a similar increase in Sersic index 
is expected for smaller, and probably more realistic, start- 
ing values of n ~ 2 for quiescent high redshift galaxies. In 
this case the maximum Sersic index might not reach as high 
values as presented here but more accurate predictions are 
beyond the scope of idealized collisionless studies like the 
one presented here. Major merger scenarios with disk-like 
progenitors do also result in an increase of the Sersic index 
but accompanied by a very weak siz e increase or even a re- 
duction in size in the presence of gas (INaab fc Truiillo 2006; 
i Robertson et alj|2009 ; ICox et alj|2006l : IHopkins et alJbood : 
Iwuvts et alj2O10h 7 Here the increase in Sersic index is partly 
driven by dissipative processes at the center of the remnants. 
Consistent with the results presented here, even considering 
subsequent 'dry' major mergers kee p the Sersic indices for 
the main stellar bod y low at n < 5. (|Naab fc Truiilldl2006l ; 
IHopkins et~ali r2009). Therefore a minor merger driven sce- 
nario - the late addition of old stars at large radii - pro- 
vides a powerful explanation for the predominance of high 
Sersic indices of pre s ent day large and massive ellipticals 
llGraham et all l200ll : iTruiillo et all |2004 iKormendv et all 



120091 : iHovod 120111 : iKormendv fc Benderl 120121 ) some of 
which are most likely the descendants of compact 'red 
and dea d' flattened high-redshift galaxies with lower Sersic 



indices (Ivan Dokkum etafl 120081 ; Ivan der Wei et al 
Bundy et alJl20ld: IWuvts et al.ll2010l: iBuitrago et al 



2009; 



2011 



Wuvtdl201]f (van der Wei et alfolllSzomoru et al.ll2012l) 

Naturally, major mergers are also expected to happen 
during the assembly of massive galaxies. There is direct 
and indirect observational evidence for this b ut the ex- 
pecte d rates are low ( ~ 0.5 — 2 since z ~ 2) (iBell et"al] 
| 006al: iMcIntosh et all |200S|; Ivan der Wei et all 120091; 
iRobaina et all 120101 : IWilliams et all l201ll : iNewman et al.l 
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120121 ; iMan et al.l |2012| ; iLopez-Saniuanl l2012h . A major 
merger will definitely dominate the late mass assembly 
history of the galaxy with a significant impact on its abun- 
dance gradients, kinematics and, eventu ally, morphology 
l|Whitelll97Sl ; iNaab et al.ll2006l ; lBoisH201ll l. However, as we 
have shown here, the expected evolution in size, surface 
density profile shapes and, eventually, dark matter fraction 
will be only moderate and not sufficient to evolve the com- 
pact high-redshift population into present day ellipticals. 
Another argument against major mergers driving the size 
evolution is the apparent absence (or very low number) 
of massive compact galaxies in the pr esent day universe 
l|Truiillo et all 1 20091 ; iTavlor et all l201Ch . Due to the low 
rates a significant number of massive galaxies will not 
have experienced any major merger since z ~ 2 and would 
remain compact (|Bezanson et al.l 120091 ). Nevertheless, a 
few massive compact have been found which eventually 
are relics of the high redshift p opulation or have formed 

recently in a si milar manner 

Jiang et all |2012| ; iTruiillo et ail 



Valentinuzzi et al, _2010: 



20121 ; iFerre-Mateu et"ail 



20121 '). Consistent with our framework these galaxies have 



low Sersic indices and our model predicts lower dark matter 
fractions compared to normal sized present day giant 
elliptical galaxies of similar mass. 

The evolution of dark matter fractions within the ob- 
servable stellar half-mass radius is also significantly differ- 
ent for major and minor mergers. In the equal-mass merg- 
ers presented here the dark matter fraction increases from 
the fiducial starting value of ~ 40 per cent to ~ 48 per 
cent for a stellar mass increase of a factor of two. This 
increase is driven by mixing processes during the violent 
merger process which lead to a real cha nge in the radial dis- 
trib ution of luminous and dark m atter jrlifz et aljl2012l . see 
also iBovlan-Kolchin et al.l l2005h . For the same increase in 
mass, minor mergers (e.g. 1:5) result in a four times stronger 
increase of the dark matter fraction (from ~ 40 per cent to 
> 70 per cent). Here the host galaxy structure is only weakly 
affected but the effective radius of the stellar distribu- 
tion increases significantly into regions which were ab-inito 
dominated by dark matter. Therefore only the ruler with 
which the galaxies are measured is changing dHopkins et all 
20091; ICenarro fe Truiillol 120091 : lLackner fc Ostrikerl I2OI0F 



Hilz et alJl20T3) . We note that the exact numbers for the 



final dark matter fractions are determined by our choice 
of the initial dark matter fraction. Our fiducial starting 
value of ~ 40 per cent might in reality be lower for mas- 
sive high-redshift galaxies but also in this case we expect 
similar trends. There is evidence for the presence of dark 
matter in present day massive ellipticals, however the ex- 
act amount is uncertain ([Gerhard et al.ll200ll ; iThomas et al.l 
120091 . l201ll : ICappellarill2012al 1. Most recent lensing observa- 
tions indicate dark matter fractions (with large uncertain- 
ties) in the range of 20-5 per cent with higher fractions 
for m ore massive galaxies (|Auger et al.ll20l"ol ; iBarnabe et al.l 
12011ft . These trends are easy to reconcile in the context of 
an evolutionary scenario where the early formation of ellipti- 
cals is dominated by a dis sipative formatio n with little dark 
matter at the center (see l|Toft et al.ll2012l ) for tentative ob- 
servational evidence) followed by a dissipationless assembly 
dominated by minor mergers naturally increasing the dark 
matter fractions. This increase would be larger for more mas- 
sive galaxies which assemble more mass by late stellar merg- 



ers (see e.g. De Lucia et al. 200£ 


; De Lucia & Blaizotl 2007; 


Guo & WhiteT 20081; Oser et al.1 


2010L 20121; lLackner et all 


2012h. However, predicted dark matter fractions depend sen- 



sitively on assumptions about the stellar initial mass func- 
tion and recent evidence for rising stellar mass-to-light ratios 
with stel lar mass leave less room for the presence of dark 
matter dvan Dokkum fc Conrovl l20ld; ICappellaril l2012b| ; 



Conrov fc van Ivan Dokkum fc Conroy||2012l : 

Ferreras et al.ll2012T ). 

More quantitative predictions, some based on sim- 
ple binary merger simulations of spheroids similar to the 
ones presented here, in a cosmological context have been 
attempted and co me to varying co nclusions; some agree 
with observations (|Ciotti et alj|2007h . some highlight pos si- 
ble tension wi t h obs e rvational trends dNipoti etaL 009bllal; 
Cimatti et ail |2012| ; iNipoti et ail l2012t iQuilis fc Truiillol 



20121 ). However, the quantitative predictive power of such 
studies is limited by construction as a realistic cosmological 
assembly of massive galaxies is very complex and many de- 
tails depend on the model assumptions about the structure 
of the galaxies and the satellites, galaxy orbits, gas fractions, 
merger rates etc. Qualitatively, however, our results are in 
good agreement with high-resolution cosmological simula- 
tions (|Feldmann et al.ll2010l ; lOser et alj|2012l ; lLackner et al.l 
l2012f ), which are supposedly the method of choice to evalu- 
ate these processes in more detail in the future. 

The predictions of this simple - dark matter as- 
sisted - minor merger driven evolution model are 
in good qualitative agreement with careful observa- 
tions of the structural redshift e v olution of ear l y-type 
gala xies llvan Dokkum et"al] l20ld: IWeinzirl etaU I2OIII 



Buitrago et al. 20 111 ; Szomoru et al.l 20121 ; McLure et al.l 
20121) . We note here that we have tested for the extreme 



case of compact progenitors. A signifi cant fraction of high 
redshift ellipticals h ave larger sizes l|Mancini et al.l |2010| ; 
ISzomoru et al.ll2012t) and would require relatively little evo- 
lution and less mergers. All the above effects are expected 
to be stronger for more massive galaxies which are expected 
to have a larger fraction of their stars acquired in accre- 
tion events. Th e most extreme cases would be centrals in 
galaxy clusters (Ostriker fc Hausmanl 1977 ; Dubinskil 1998 ; 



De Lucia fc Blaizotl l2007i ; iRuszkowski fc Springell I2OOS 
Laporte et al. 20121 ) 



If we are right, then certain definite and testable pre- 
dictions can be made about the properties of the outer parts 
(r > r e ) of normal giant ellipticals: Since they are the result 
of accretion of low mass systems they should have consid- 
erably lower metallicities than the inner parts. Steeper gra- 
dients at large radii then mi ght imply a f ormation history 
dominated by minor mergers jWhitdll978l ; lKobavashill2004l , 
see however iPipino et alj|2010] for an alternative view) . The 
stellar orbits at large radii should become radially biased 
i|Hilz et alj|2012h and the dark matter fractions of massive 
ellipticals should increase rapidly at large radii. 
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